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We study the von Neumann and Renyi bipartite entanglement entropies in the thermodynamic 
limit of many-body quantum states with spin-s sites, that possess full symmetry under exchange 
of sites. It turns out that there is essentially a one-to-one correspondence between such ther- 
modynamic states and probability measures on CP 2s . Let a measure be supported on a set 
of possibly fractal real dimension d with respect to the Study- Fubini metric of CP 2s . Let m 
be the number of sites in a subsystem of the bipartition. We give evidence that in the limit 
m — > oo, the entanglement entropy diverges like | logm. Further, if the measure is supported on 
a submanifold of CP 2s and can be described by a density / with respect to the metric induced 
by the Study-Fubini metric, we give evidence that the correction term is simply related to the 
entropy associated to /: the geometric entropy of geometric quantum mechanics. This extends 
results obtained by the authors in a recent letter where the spin-g case was considered. Here 
we provide more examples as well as detailed accounts of the ideas and computations leading 
to these general results. For special choices of the state in the spin-s situation, we recover the 
scaling behaviour previously observed by Popkov et al., showing that their result is but a special 
case of a more general scaling law. 
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1 Introduction 



The study of the entanglement entropy of extended quantum systems such as quantum spin 
chains has attracted much attention in recent years (see e.g. the recent special issue [1]). The 
entanglement entropy is an interesting quantity for many reasons. Among other features, it gives 
information about the quantum state of a system, particularly about the amount of entanglement 
that it may store, and it exhibits "universal" behaviour, for example in systems at or near 
conformal critical points. 

There are different types of entanglement entropy but a widely studied quantity is the von 
Neumann entropy. Formally, it may be defined as follows: consider a quantum system which 
we partition into two susets A and A and define its Hilbert space H = %a ® "Ha- Suppose the 
system is in a pure state | ip) . The bipartite entanglement entropy S is the von Neumann entropy 
associated to the reduced density matrix p A of the subsystem A, 

P A = ^n A mm, (i.i) 

given by 

S = -Tr HA (p A logp A ). (1.2) 
A generalization of the von Neumann entropy is the Renyi entropy which is defined as 

S. = '° g(Tr "^», (1.3) 
1 — n 

and depends on an additional parameter n. In the limit n — > 1 it equals the von Neumann 
entropy. The quantity p A may be interpreted as the reduced density matrix of subsystem A in 
a 'replica' theory consisting of n copies of the original model. Hence, one often refers to n as a 
'number of copies'. 

For quantum systems at conformal criticality, it is well know that the von-Neumann and 
Renyi entropies of a large block of spins of size to diverge as 

c , c 



log to and S n ~ - (l + n log to, (1-4) 
3 6 

where c is the central charge of the conformal field theory which describes the quantum critical 



point [21 |3l 0]. The first equation in (1.4) is arguably the best known result in this field of 
research and one that has been successfully tested for a great variety of theories, both through 
numerical and analytical methods. 

Less known is the fact that the entanglement entropy can also reveal very interesting in- 
formation about the state of quantum systems away from conformal criticality. An important 
body of work exists where the entropy of quantum spin chains with random interactions has 



been studied (see e.g. the review [5]) leading to scaling laws for the entropy similar to (1.4) but 
where the coefficient ^ is often replaced by an irrational number. Also, a number of relatively 
recent works [6l 171 IH1 l9l 110} [TTj involving for instance non-critical quantum spin chains (i.e. spin 
chains whose thermodynamic limit is not described by conformal field theory) or models with 
long-range interactions, have found a variety of behaviours for the entropy of large blocks of 
spins. Our own work [11] has, for the first time, provided a geometric interpretation as well as a 
generalization of some of those results for the spin-^ case. In [TT] we studied the entanglement 
entropy of an infinite spin-^ chain whose ground state is an infinite linear combination of basic 
permutation symmetric, zero entropy states. We found that both the von Neumann and Renyi 
entropies of a large block of spins diverge logarithmically with the size of the block to as 

S,S n ~ ^logTO + 0(1) with d€[0,2]. (1.5) 
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In particular, and in contrast to the results in one-dimensional conformal critical systems, we 
obtain the same behaviour for the von Neumann and Renyi entropies, as well as the single-copy 
entropy \12\ I13j . that is, the limit linin^oo S n (we will see that the 0(1) term is finite in this 
limit). 

A particular choice of the state (known as a Dicke state [13]) leading to d = 1 is the 
ground state of the ferromagnetic XXX chain, and was studied in [6J, [10] . Also, in [8], this and 
other permutation symmetric states (certain linear combinations of Dicke states) were studied 
as ground states of the Lipkin-Meshkov-Glick (LMG) model |15|. [16] I17Q In |llj we argued that 
the quantity d can be any real number in the given interval and that it represents a dimension 
which characterizes the geometry (fractal or otherwise) of the quantum state. For spin-^ all 
allowed quantum states of the system have support on the Bloch sphere of geometric dimension 
d = 2. In other words, all states of the system can be expressed as linear combinations of states 
|Vv) labeled by a vector v on the Bloch sphere and as the length of the chain tends to infinite, 
every vector in the Bloch sphere corresponds to a basis state of the system. In particular, we 
may consider a quantum state associated to the 'embedding' of a fractal object in the Bloch 
sphere, such as the Cantor set shown in Fig. 1 of Section 6.1 and recover the fractal dimension 
of such object by computing the bipartite entanglement entropy. 

In the current work we put the results above on firmer mathematical footing as well as 
showing that the method employed in [11] can been successfully generalized to higher spins s 
and general permutation symmetric states. We identify the geometry of the support of the 
states of such a system and find that the Bloch sphere is naturally generalised to the projective 



complex space CP 2s . We show that for such states the entanglement entropy behaves as in (|1.5|) 



with d G [0,4s], and that the explicit form of the correction term 0(1) in ( |1.5[ ), which is given 
in the next section, naturally singles out the Study- Fubini metric of CP 2s . The present work 
provides mathematically precise statements (see below Theorems [I] and [2j and Conjecture [3]), 
as well precise and compelling derivations (if not rigorous proofs) of them. 

The paper is organized as follows: Due to the technical nature of some of our results we have 
decided to start the paper by introducing the notation and main terminology we will be using 
in the remainder of the work. We then summarise, discus and set our main findings in context, 
comparing them when possible to previous existing work. We do all of this in Section [2} In 
Section [3] we provide a proof of one of our main results (Theorem 2 stated in Section [2]). This 
result is a general formula for the Renyi entropy based on the use of cyclic permutation operators 
introduced in previous work by the authors. In Section [4] we consider the Thermodynamic 
limit of permutation symmetric states, that is we find a suitable representation for permutation 
symmetric states when the length of the spin chain tends to infinity. This representation is 
dictated by the need to ensure normalization and orthonormality of permutation symmetric 
states in the thermodynamic limit. In Section [5] we compute the bipartite entanglement entropy 
of certain linear combinations (both finite and infinite) of permutation symmetric states in the 
thermodynamic limit. In each case, we extract the leading behaviour of the entropy as the size of 



the block tends to infinity. In Section 
entanglement entropy as in Section 



we first perform a similar computation of the bipartite 
but we specialize to s = | and consider a very special 
linear combination of permutation symmetric states whose support is the Cantor set projected 
onto half of a great circle on the Bloch sphere. We then provide a general argument on how the 
entropy of large blocks may scale for linear combinations whose supports have fractal geometries. 
Finally, we close the paper with some conclusions and outlook on Section 7. 



*Note that it was found in [TS] that in the LMG model the same critical behaviour is observed for the von 
Neumann and single-copy entropies. 
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2 Main results and discussion 



2.1 Main results 

Consider a Hilbert space composed of N copies of spin-s representation spaces C D+1 , D = 2s: 

H N =(C D+1 f cN 

Let Vn be the subspace of Hn formed by all vectors that are symmetric under all permutations 
of sites. It is a simple matter to see that one can describe basis elements for Vn by specifying the 
numbers Nj of sites with ^-eigenvalue equal to s — j, j = 0, 1, . . . , D, under the sole condition 
Y2j=o Nj = N. Any such basis vector for Vn is a suitably normalized linear combination with 
equal coefficients of all vectors with fixed Nj: 



^(iVciVi,...,^)) 



N lNi\...N D \ 



^ |cr^ Vo- - -VQVi • • •Vj ... Vp . ..V/j 



(2.1) 



where the sum is over permutations of the vectors Vj, with i = 0,1, ... D. The vectors Vj are 
complex D + 1-dimensional column vectors with one single non- vanishing entry of value 1 at line 



i (starting from line 0). For every fixed choice of the parameters No, • • • , Njj the state (2.1 ) is 
an eigenstate of the total spin operator S z = Y2i=i a \ with eigenvalue \x = YliLo( s ~ We 



will refer to the basis vectors (2.1 ) as elementary vectors (these are also called Dicke states 

Consider the particularly simple permutation symmetric vectors where all sites are in the 
same local state described by v: 

K) = |v,v,...,v). (2.2) 

It turns out that it is also possible to write any permutation symmetric vector as a linear 
combination of IVv) 8 - These are more useful for our study of the entanglement entropy, because 
the vectors (2.2) are not entangled, contrary to the elementary vectors (2.1). They are the 



only permutation symmetric vectors that are not entangled. We will refer to (2.2) as zero- 
entanglement vectors. 

In order to see that the zero-entanglement vectors span Vjy, let us use the following parametriza- 
tion for normalized complex column vectors, 



V P, a :-- 



( a e ie ° \ 
aie iei 

\ a D e ie » J 



D 



with 



Oj = 1 and a,j € K Vj. 

j=0 



(2.3) 



We use 6 
for i = 0, 



(8q, . . . , Od) with 9i £ [0, 2ir] for i = 0, . . . , D and a := (ap, . . . ,an) with G [0, 1] 
, D. It is a simple matter to see that 



v(9,a)l 



N 

E 

N ,...,N D =Q 



D 

n 



a 3 



6 N,J2f =0 N J 



|*(JV ,JVi,...,jyb)). 



(2.4) 



The inversion of (2.4) is obtained by performing the following operations: first, we multiply both 



sides of the equation by exp(— i X^q 1 finally we integrate both sides of the equation over 
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the variables 9q, 



, Qd—\-> normalizing each integral by a factor of 1/2-7T. We obtain: 



1 



2tt 



(2vr)^ ,„ 

2- 



d D 9e 



v(0,a) 

d D ^Ejio 1 e J {N J -N ) e ie D N D8 



(27T) 



N 



N ,...,N D =0 

v ^r, v .|*(JV ,iVi,...,JVb)) 




where we employed the shortcut notation: 



2tt /.27T 

d D ej{6) : = / 
JO 



2tt />27r 

^1.../ de D - l f(6 Q ,...,e D „ 1 ). 

JO 



(2.5) 



(2.6) 



Employing <5 a .b = ^ J 2?r e^ a we can carry out the integrals on the r.h.s. to obtain: 

i / ° „/m\ i 



\^(N ,N U ...,N D )) 



D 

n 



' N! \iS 



(2vr) D 



2/r 



d D 9e i=° 







(2.7) 



We have integrals rather than strictly a linear combination, but it is obvious that these integrals 
are well defined (e.g. as integrals on C N ( D+1 ^) and give a vector in Vn- It is interesting to note 



that (2.7) is in fact independent of the choice of the parameters a. These can in principle be 



fixed to any particular values, as long as the normalization condition (2.3) is met. Some choices 



however, are more appropriate for certain situations; an interpretation for these parameters will 



become apparent later, when the thermodynamic limit is considered. From (2.7) it is easy to 



show that the state |^(iVo, N\,..., Ne>)) still has norm one as it should, using the overlap 

N 

<^v ( y VV(£,,)> - I y \ <'i~i<'" H '~ H ' ' I , (2.S 



vi=o 



which is immediate from the definitions (2.2) and (2.3). 



Equation (2.7) means that the zero-entanglement vectors (2.2) span Vn- Clearly, from (2.7) 



we may as well fix both the length and the overall phase of v in (2.2), and the vectors will still 
span Vn- Hence, we may see v as homogeneous coordinates for CP , and we find that vectors 
parametrized by points in CP D span Vn- This is an over-determination: as we mentioned, it 
is possible to fix a, and it would in fact be possible to take only a finite number of points in 
CP 15 in order to span Vn, although we will not need this here. Note that CP 2s is the geometric 
space of quantum spin-s states, studied from this viewpoint in the context of geometric quantum 
mechanics |19| , I2Q [, [2Tj . We will discuss this connection below. 

Let us now consider the limit N — > oo of Vn- This limit of course must be taken with 
care. In particular, we need an appropriate topology. Taking the somewhat complicated limit of 
vectors in the Hilbert space is not particularly useful for our purposes. Instead, we will consider 
quantum states seen as linear maps on End(%jv) : 



End(ftjv) 
O 



101 



(2.9) 



and take limits on these linear maps. We will denote by Pat the set of quantum states (linear 
maps) corresponding to vectors in Vn- That is, the set Pyv is the set of pure permutation 
symmetric quantum states on iV sites. 
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For quantum states corresponding to vectors in T-Ln in general, instead of Vn, one would 
require, in order to define the large- iV limit, an embedding of Hn into Hn+i for every N: a 
prescription as to where the sites are being added (in more mathematical terms, in order to 
define the limit set of operators, one would have to construct a direct system on the algebras 
of linear operators on Hn, and take the direct, or inductive, limit). However, since we are 
looking at permutation symmetric states, there is a canonical embedding. This goes as follows. 
Any operator O 6 End(%Ar) is a linear combination of products of operators A = A\ ■ ■ ■ An 
factorized on the tensor factors of %n- By permutation symmetry, when evaluating if (A) with 
if) € Pjv") the factors acting nontrivially (i.e. different from the identity) can always be gathered 
contingently on the sites, say, at positions 1,2,3, . . ., and in any order. Hence, for a factorized 
operator A, the only information that we need in order to evaluate if (A) is the set of nontrivial 
factors, as matrices on C D+1 (we do not need the information of the site on which they act). The 
same information can be provided to evaluate if' (A) for any \if') £ Vmi as long as the number 
of nontrivial factors is less than or equal to M. This, along with linearity, gives the canonical 
embedding if M > N. 

Let us refer to any linear operator O 6 End(%Ar), for any finite N = 1, 2, 3, . . ., as finitely- 
supported (following the terminology of [10J ) . Let if N £ P^r, N = 1,2,3,... be an infinite 
sequence of quantum states. We will say that lim7v-s>oo iPn exists in the local- operator topology 
if limTv-s-oo ^jv(C) exists for every finitely-supported O (with the embedding described above). 
This defines a linear map on the space of finitely-supported operators, if) = limjv->oo iPn with 
if{0) = limj^^oo ij) N (0). Let us denote by P the space of all such linear maps, that occur as 
limits in the local-operator topology. We may see this as the space of all permutation symmetric 
quantum states in the thermodynamic limitj^] As described above, if if) E P and A is a finitely 
supported operator that is factorized, then if (A) only depends on the non-trivial factors of A. 
Also, if) has the properties that if(l) = 1 and that if)(0)* = if)(0'), as usual for quantum states. 
Note that an explicit metric inducing the local-operator topology (in a slightly more general 
context) was given in |10| . 

For any fixed v, the N — > oo limit of the sequence of quantum states corresponding to zero- 



entanglement vectors \if) v ) defined in (2.2) is obviously convergent in the local-operator topology. 
By a slight abuse of notation, we will denote it by if v £ P, and refer to it as a zero-entanglement 
state. Our first result is as follows. 

Theorem 1 Let if £ P. Then there exists a probability measure /u on CP D , 

d//(v) = 1, 



such that 

if)(0)= [ d^)if v (0) (2.10) 
for every finitely- supported O. Further, every probability measure gives rise to a unique if G P. 

Note that we expect the measure to be unique, since in particular CP 13 is compact (there is 
obviously unicity up to the weak equivalence under equality of all averages of finitely-supported 
operators). 

Basically, what is happening is that all of CP D becomes necessary to span Vn in the large- N 
limit. In fact, CP D becomes a set of basis elements, and further (if v \0\if v i) — > whenever v 



^The set of limits limjv^oo V'jvi with ij) N sequences of pure permutation symmetric quantum states, is expected 
to give rise, from the viewpoint of local operators, to all thermodynamic permutation symmetric quantum states, 
pure or mixed; but we do not have a proof of this. 
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and v' are not colinear, for every finitely-supported operator O. Hence, any linear combination 
I^n) = XX^lVv) gives rise, in the large-iV limit, to the quantum state XX l c v| 2 Vv formed 
purely out of ip v , because the cross-terms vanish. The sum may contain infinitely many terms, 



and the more precise way of describing this is using a measure: XX 



i — ^ 



f CpD d/i(v). 



In the present paper we will not present a complete proof of Theorem [TJ but in Section [4] we 
will give various derivations which provide a good intuition for (2.10). The complete proof of 
Theorem [T] will be presented in a forthcoming work. 

The quantity of interest to us is the entanglement entropy (1.3). Our next result shows that 
the above formalism is useful to evaluate the entanglement entropy in the thermodynamic limit. 



Theorem 2 Let ip N S Ptv ; N = 1,2,3,... be a sequence of quantum states converging as 
N — > oo to ip 6 P in the local operator topology. Then for every integer n > 1, the Renyi entropy 
of the states ip N , with a fixed number m of sites in a subsystem of the bipartition, converges as 
N -+ oo to 



Sn 



1 



1 



n 



log 



(Cp-D)xn 



n dfi(v a ) n 



.1 



(2.11) 



^a=l / \o;=l 

where fi is the measure of Theorem^ andvX is, as usual, the transpose of the complex conjugate 
of the column vector v a . Here we use v n+ i := vi. 



Note that on the right-hand side of (2.11) the integrand indeed lies on (CP D ) xn ; in particular, 
the phases of v Q are unimportant thanks to cyclicity of the product. We provide a proof of 
Theorem [2] in Section [3] 

There is a canonical metric on CP D : one considers the metric on the sphere S 2D+1 = 
{(z , zi,...,z D )e C D+1 : |z | 2 + • • • + \z D \ 2 = 1}, and takes for CP D = S 2D+1 /S 1 the quotient 
metric induced by simultaneous phase shifts on all coordinates z,. This is the Study- Fubini 
metric (22^ [23], [2^1 [25] (also referred to as the Bures-Uhlmann metric in the context of geometric 
quantum mechanics |26[ [2"T] ) , giving in particular a differentiable manifold structure to CP D of 
real dimension 2D. Our next results, which are the main results of this paper, show that the 
entanglement entropy is naturally associated to this manifold structure and in particular to this 
metric. 



Conjecture 3 

I. Let the measure fi of Theorem^ be supported on a subset U of CY D with (possibly fractal) 
dimension d with respect to the Study-Fubini ( Lip shitz-) manifold structure. Then the large-m 
behaviour of S n is given by 

S n ~ ^logm + 0(l) (2.12) 

for every real n > 1. Note that < d < 2D = 4s, and all cases are possible. 

II. Assume that the subset U where \x is supported is a submanifold ofCP D . Assume further that 
\x is absolutely continuous with respect to the volume element on U induced by the Study-Fubini 
metric on CP D . That is, there exists a density /(v) such that d/j,(v) = <i d v/(v) where d d v is 
the volume element on U . ThenH 

S n ~^log^ + r ^log(Vi y>v/(v)") +o(l). (2.13) 

''Note that for spin-| this result was first given in [11], albeit containing some typos. 
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In Part I of this conjecture, we expect only the Lipshitz-manifold structure to be important, as 
the fractal dimension is invariant under bi-Lipshtiz transformations. 

Results (2.12) and (2.13) are derived from (2.11) taking n > 1 integer. However, they make 
sense as well for real n > 1. We conjecture that they are the correct asymptotic results for the 
Renyi entropy for all real n > 1. A general calculation as well as various examples for Part II of 
Conjecture [3] are provided in Section [5} and for Part I, in Section [6j 

With this conjecture, we find in particular for the von Neumann entanglement entropy, 



Si 



d em 
- log — 
2 B 8vr 



J d d v/( V )log/(v) + o(l) 



(2.14) 



We observe that the second term in this expression is the so-called geometric entropy of geometric 
quantum mechanics. This is the entropy associated to the probability density /(v). 

Another interesting limit is the limit n — > oo, which gives the single-copy entropy Soo = 
— log Ai where Ai is the largest eigenvalue of the reduced density matrix. A simple saddle-point 
analysis gives, in the case where / max := max v /(v) < oo, 



Sr, 



d m 
- log — 
2 6 8vr 



log(/max) + o{l). 



(2.15) 



Note that in all cases, the leading large-m behaviour is the same, in contrast to the results 
for critical one-dimensional quantum chains. 

Note that in (2.13) and (2.14), the Study- Fubini metric is singled out: it is the density 
/(v) with respect to the Study- Fubini-induced volume element that occurs, and a different 
choice of density would lead to different-looking, more complicated formulae, since /(v) n and 
/(v)log/(v) are not densities (in particular, there would be n-dependent extra factors under 
the integral). 

To be more precise, the normalization of the metric that we choose is such that, in particular, 
the sphere CP 1 has radius 1. The distance function is, in homogeneous coordinates |21| I28j. 



D(v,w) = 2arccos 



V I W W I V 



V I V W I w 



(2.16) 



We can write down the Study-Fubini metric and volume element explicitly. Consider the coor- 
dinates 6i 6 [0, 27r), i = 0, . . . , D — 1 as well as the spherical coordinates i = 0, . . . , D — 1 on 



the sphere {(do, ai, . . . , ao) £ M. D+1 : 



D 



1}, 



a 



COS ( 



'0) 



COS i 



n 

k=0 



D-l 



8111 ( 



for 



1, . . . , D — 1 and ajj = J^J 



sm< 



k=0 



The at are nonnegative, hence 



< 



< 



7T 



for 



D-l. 



In terms of this parametrization, the Study-Fubini metric is 



fc-i 



4* fcJ - n 



sin 



(2.17) 



(2.18) 



(2.19) 



leading in particular to the volume element 

d 2D v = 2 D J] sin 2,; 

p=0 



D-l k-1 

nn 

fc=0 p=0 



sin 



p=0 



d D 0d D c 



(2.20) 
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2.2 Discussion 



Geometric quantum mechanics. Our results have connections to geometric quantum me- 
chanics. Indeed, in that context, one also uses probability measures fi on CP D , in order to 
describe the mixed states of statistical quantum mechanics. As is well known, this description 
is rather redundant from the viewpoint of linear quantum mechanics: one can always diagonal- 
ize the density matrix, hence choose a measure supported on the points corresponding to the 
eigenvectors. An important difference, in our set-up, is that there is no redundancy: different 
measures \x correspond to different thermodynamic permutation symmetric quantum states. 

An natural quantity in the context of geometric quantum mechanics is the geometric entropy 
— f d d v /(v) log /(v) associated to a Study-Fubini density /(v). This quantity, in particular, 
partially lifts the redundancy mentioned above, hence its meaning, in the context of linear 
quantum mechanics, is not always clear. Further, one can ask about the maximization of the 
geometric entropy under the condition of a fixed averaged energy, controlled by the temperature 
T. It turns out that it does not lead to the usual thermal density matrix of statistical mechanics, 
where eigenstates of the energy appear with Maxwell-Boltzmann probabilities. Rather |21j . it 
leads to the density 

e -H(y)/T 

WV) = J cpP ^v e -»W/T - < 2 ' 21 > 

where -ff(v) is the average of the energy in the quantum state pointed by the homogeneous 
coordinates v of CP D . It can be shown that this is not simply a redundant expression for 
the standard thermal density matrix; it gives a different linear mixed state [21 1. Our set-up 
provides an appealing physical interpretation for /*. It is the distribution for thermodynamic 
quantum states obtained by maximizing the entanglement entropy of a large subsystem, under 
the condition of a fixed average energy, and under the extra condition that the system be 
restricted to the subset of permutation symmetric states. 

Creating entanglement. A standard way of creating entanglement in the context of quantum 
information theory is by symmetrizing a vector formed by factorized qubits. The resulting 



states are usually referred to as Dicke states [14J. Our results (2.12), (2.13) provide an upper 
bound for the entanglement that can be stored by this method in a large subsystem of m spin-s 
qubits singled out from a much larger system of spin-s qubits. The entanglement can grow 
logarithmically (as already found by Popkov et al. [7]), with a coefficient that can be as large as 
D = 2s, but not larger. Further, when this bound is achieved, the entanglement is maximized 
by maximizing the geometric entropy associated to the coefficients obtained by expanding in 



the zero-entanglement vectors (2.2). 

Note that in [7J it was stated incorrectly that the logarithmic growth for the entanglement 
of pure permutation symmetric states couldn't have a coefficient larger than s. It was however 
remarked that for mixed states, the coefficient could be as large as 2s. Here we have not analyzed 
mixed states, but we expect that the set P, although obtained as thermodynamic limits of 
pure permutation symmetric states only, contains also all mixed permutation symmetric states; 
whence our results being consistent with Popkov et al.'s remark. This is based on the known fact 
that any mixed state on N sites can be reproduced using a pure state on 2N sites. Hence in the 
large- N limit, from the viewpoint of finitely-supported operators, there is no distinction between 
mixed and pure states. This argument is much more subtle when one restricts to permutation 
symmetric states, because the doubling procedure of a mixed permutation symmetric state on 
N sites does not, in general, produce a pure permutation symmetric state on 2N sites. It would 
be interesting to analyze this situation in more depth. 
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Comparison with previous results. In [7J, the limit N — > oo was taken on the entan- 
glement entropy (in fact, on the reduced density matrix) associated to the elementary vectors 
|\P(iVo, N%, . . . , No)) by keeping fixed ratios Nj/N = pj (more precisely, by "keeping fixed ratios 
Nj/N = p^ 1 we will mean taking Nj = Nj(N) such that linijv^oo Nj(N)/N = pj). Clearly, we 
have X^jLo^i = We wm show that this limit exists in the local operator topology on the 
elementary quantum states *$>(Nq, Ni, . . . , Nu) themselves: 

y(p , Pl ,...,p D ) := lim *(JVo,JVi,...,iVb). (2.22) 

iV-»oo 

We will show that in this case, the measure \i above is supported on the submanifold specified 
by a% = y/pl, i = 0, 1, . . . , D, and that on this submanifold, we have 

dfx(v(9,a)) = —-^d9 ---dO D -i, 9 D = 0. (2.23) 

We provide the derivations (but no mathematically rigorous proofs) in Subsection |4.2| 



This provides the interpretation of the coefficients a% that were arbitrary in (2.7): they must 



be fixed to the square- roots of the ratios Nj/N in order for the thermodynamic limit to make 



sense on (2.7). In particular, the large-iV limit considered in is not the most general one. 
The leading logarithmic divergence found in [7] was slogm, because the chosen submanifold is 
of dimension d = D = 2s, half of the real dimension of the CP D manifold. 

For more precision, let us note that the metric induced on the submanifold ai = yfpi is gej,e k 
(2.19). It has volume element 

d D v = 2 D Y[ ai d D 9. 

i=0 

Hence, we see that the conditions of Conjecture [3}ll are satisfied, with 

I 



(4^n?= ^ 



Therefore, equation (2.13) gives us, in terms of s = D/2 and p% = cq, 

slo n 1 ° 

S n ~ slog2vrm + + - logfbi + o(l). (2.24) 

n — 1 2 

i=0 

For n = 1 we recover the result obtained in [7J, namely the large m behaviour of the von 
Neumann entropy: 

1 ° 

5i~slog27rem + -logJJpi + o(l). (2.25) 

i=0 

Finally, let us note that for spin 1/2, it was found in [8] that the ground state of the Lipkin- 
Meshkov-Glick model, for certain values of the parameters where it is not a pure Dicke state, 



has an entanglement entropy with the logarithmic behaviour (1.5) with d = 2/3. It would be 



very interesting to understand this state from a geometric point of view. 

3 Proof of Theorem [2] 

A proof of this result for s = \ was sketched in the letter [TT] . Here we will present an extended 
and more general proof. The starting point is the expression of the entropy in terms of local 
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cyclic replica permutation operators 71 which act on site i of a quantum spin chain and cyclicly 
permute the spins of n replicas of the model at that same site i. These operators were introduced 
in |10j where we also showed that the Renyi entropy of a block of (not necessarily consecutive) 
spins of size m is given by 



1 — n 



log 



where ip® n is the n-th tensor power of the state i/j and 

T A = Y[Ti and \A\ 



m. 



(3.1) 



(3.2) 



Employing the language of functionals as in the previous section, we can express the correlation 
function above simply as 

„_*(*-C7i))_ 



From |10j we also know that 



1 — n 

n D+l 

% = Tr aux [ ]J E E ^ E Z, 
,a=l ei,£2=l 



(3.4) 



where E eie2 are the (D + l) x (D + l) elementary matrices, where all elements are zero except for 
a 1 in row ei, column €2, and the indices aux, a and i refer to an auxiliary space V aux = C D+1 , 
the copy number and the site, respectively. Since 7a, for every finite set of sites A, is a finitely- 
supported operator, and since the entanglement entropy is expressed as the evaluation of this 
operator in a quantum state, then by Theorem [T] the thermodynamic limit of the entanglement 
entropy exists and can be expressed using ( 2.10[ ). 

We can now write the Renyi entropy of the infinite chain by recalling (2.10) and applying it 
to the replica case and the operator Ta 



(CP D )> 



n^(v a ) (g)w Q (TA) 



\a=l 



a=l 



with 

n 

(8)^v Q (TA) = n Tr au: 

a=l i£A 

where we used, in an obvious notation, 



n e E ™^«A E %i) 1 > 

,o=l ei,€2=l 



The quantity ip Va ^(E^f tl ) is independent of the site i, and from (2.2) and (3.7), 

i> Va ,i( E %k)=<ei V <x,ei- 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



Substituting into (|3.6|) and tracing over the auxiliary space we find 
Tr 

aux 



n D+l 

n e k 

.0=1 ex,£2=l 



aux * 

62 U a,€2 Ua > e i 



D+l 

E ^1^1,^2,^2,, 
ei,...,e n =l 
n 

Yl v l ■ v «+i • 
0=1 



v n,e n V nei 



(3.9) 
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Therefore 



(CP D )> 



n ^( v «) n v « ■ vq+i 



(3.10) 



\a=l 



\a=l 



which when substituted in (|3.3|) gives the result (|2.11l) for vj, • v Q = 1. 



4 Thermodynamic limit of permutation symmetric states 

The complete proof of Theorem[T]is rather involved, and slightly beyond the scope of the present 
paper. We instead present, here, a proof of (2.10) in the case of finite linear combinations of 



vectors m fixed directions v, and a derivation (not of mathematical rigor) of the expression 
( |2. 10 ) , (2.23) for the states (2.22). The latter are the thermodynamic limit N — > oo of the states 
fy(No,Ni, . . . ,Nd) associated to (2.1), with the condition that the ratios Nj/N = pj be fixed. 



This not only will give us some intuition as to the way the measure \i of Theorem [T] appears, 
but also will be immediately relevant to the analysis, within our framework, of the results of 
Popov et al. [ZJ. We will also consider the cases of finite linear combinations of the vectors 
|^(iVo, N±, . . . , Ne>)), held in the thermodynamic limit at different sets of ratios {pj}- 



4.1 Finite linear combinations 

Let V be a finite set of complex vectors v of unit length v T v = 1, such that no two of them are 
colinear (i.e. related to each other by an overall phase). Let us consider |V>at) = J2 v <=v °v,N |VO £ 
T-Ln-, with the condition that limTv^oo c v ,at = c v exists. For every finitely-supported operator O, 
which acts nontrivially on I sites, we have, for all N > I, 



<W|0|Vw> = (^vl^w>Xv,w(0) 



where the factor 



Xv,w(0) 



(4.1) 



(4.2) 



is independent of N. In particular, we note that ^ vw (0) is a polynomial in the components of 
v and w (hence continuous), and that 



Xv JO) = w(0) 



(4.3) 



(recall that by our convention, ^ v is the limit iV — > oo of the quantum states associated with 
the finite-iV vectors IVv))- Clearly, we have (V'vlV'w) = (v'w)^. The maximal value of |vTw| 
occurs when v and w are colinear, hence equal to each other if both are in V. This maximal 
value is 1. Hence, 

1 (v = w) 



lim (V'vlV'w/ -in 
JV-S-oo I (v / Wj 



(4.4) 



whenever v, w 6 V. This yields 



lim <V;vP|Viv) 

N^oo 



lim V c; iV c w^ 7v(V' v |V'w>Xv.w(C , ) 

N^oo * — ' 

v,w£ V 

E l C vl^v(0) 



(4.5) 



vev 
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which shows (2.10) in the cases of finite linear combinations, with 

f ^(v)F(v) = ^|c v | 2 F(v). 



(4.6) 



vev 



That is, the measure is supported on the finite set V, with weights |c v | 2 : 

lim tp N = V IcvlVv- 



(4.7) 



4.2 Elementary states 



We now consider the thermodynamic limit of the elementary states ^(Nq, N\, . . . , Nd) associated 



to th e elementary vectors (2.1). In order to apply a similar methodology, we use the expression 
(2.7). As we said, the ratios ^ - : ' 1 : v 



N 



Pi remain finite, and since 

N = Ei=o N i [t follows that 
Yli=oPi = 1- For simplicity of the derivation, we will assume that pj ^ for all j, but this is 
not essential and the final result holds for general pj. 

The intuition behind (2.10) is that, as is clear form the previous subsection, the vectors 
|Vv) become orthogonal, at different values of v, in the limit N — > oo. Since in this limit 
all directions of v are independent, it is natural to expect that in regular enough situations, 
like that of ^(Nq, N\, . . . , Nd), we should normalize the resulting vectors with a delta- function 
normalization. If, for instance, only the angles 9j may vary, then we will have vectors |^ , v (9,a)) 
representing the thermodynamic limit of the vectors \ip v (e, a )) taken in such a way as to ensure 
that they form a basis of factorizable, orthonormal states under the following delta-function 
normalization 

D-X 

(*vd,«)l*v(i ) a))=II 27r ^i-^)- ( 48 ) 
3=0 

The vector \^ v (e.a)) i s > formally, a limit of the form lirn/v->oo Rv(N)\tp v ) for an appropriate 
renormalization factor R V (N). Note that ( |4.8| ) implies, for every finitely supported operator O, 

<*v(^l°l*v(0,a)> 



D-l 
3=0 



(4.9) 



(O) (see (4.2)) is independent of N, a fact which we used 



Here, the factor Xy_{°) := Xv&a),v(&a) 
on the third line in order to evaluate the limit. On the last line, we used (|4.3l). 



It turns out that in the thermodynamic limit we can express (formally) the vector 
|*(po,Pi,...,Pu)> := ArAr lim \*{N ,Ni,...,N D )) 



as 



|*(p ,Pi, • ■ • ,Pd)) 



(2tt) 



D 



2tt 



d D 9\V 



v(9,a), 



(4.10) 



(4.11) 
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where the integration is on 9q, . . . ,0d-i with Op = fixed, and the aj are fixed in terms of the 
Pj as 

a 3 = VpJ- ( 4 - 12 ) 



This implies the expression (2.10), with integration on the submanifold given by the fixed aj 
and measure d D 9/(2ir) D as claimed in ( |2.23| . Indeed, 

9(po,Pi,...,p D )(0) = (*(po,pi,...,pd)\0\V(po,p 1 ,...,pd)) 



(2*) 



2D 



277 



d D ei O<* v(8fB) |0|*v(aa)) 



2?r 



(2^r^^ (o) 



(4.13) 



where we used (4.9). 

We now present a derivation to find the coefficient R V (N) such that (4.8) holds, and to show 
that CTil holds with Cl2|). 



Let us start by considering the numerical coefficient in (|2.7|) 

D 




Nj\ 



D 



(V2^N) D l[p^ 

3=0 



/Pj 



(4.14) 



where we have used Stirling's formula A\ ~ v / 2vrA^ r for A large. We can deduce from (2.7), 
(2.8) and (4.14) what the state |^ , v (e,a)) must be. Since \^(po,pi, ■ ■ ■ ,Pd)} must have norm one, 
it follows that 



N 



D D 



(27T) 



D 



3=0 



2Nj r 2n 



d D 9 



2tt 



D 



N 



, (4.15) 



which may be rewritten as 



2tt 



d D e x 



(2k) 



D 



D 



3=0 



.fc=0 v ° fc 



D 

£ 



A? 



. (4.16) 



This can be written as the integration of a multi-dimensional delta-function. Indeed, if the sum 
under the iV-power on the r.h.s. has modulus greater than or equal to one for 9 = 9 and less 
than one for 9^ 9, then in the large- N limit the integrand will be supported on 9 = 9, hence 
will be a multi-dimensional delta function. This condition is only met if we set aj = ^fp- for 

all j = 0, ... , D; then for 9 = 9 the terms in the sum add up to one because Y2f=oPj = 1, and 
otherwise they add up to a number whose modulus is less than one. This implies that 



(2itN) 



D 



D 



3=0 



D 



N 



J=o 



iV->oo 



D-l 



Y[ 2ir5(9j - 9j). 
3=0 



(4.17) 



Thus, the vectors \^ v ($ a )) are °f the same type as their finite N counterparts with 



D 



|* V(M > = Jim (V2^N) D H ^a]e 

TV— >-oo 

3=0 



'#v(fl,a)>> 



(4.18) 



that is, Rv(e,a) (N) = ( 2 ^ N ) s Uf=oVaj e ^ 3Nj - From (j 2 - 7 ^ this implies (|4.1l|) and (|4.12|). 



,-JV, 
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4.3 Finite linear combinations of elementary vectors 

We consider vectors 



(4.19) 



where each element in the set Qn of multiplets N = (Nq,N\, . . . ,Nd) is associated to a N- 
independent multiplet p = (po,Pi, ■ ■ ■ ,Pj) € Q with YljPj = 1> such that limjv^oo Nj/N = pj. 
We assume that liniTv^oo cat = c p exists. 

The limit N — > oo of the state ip N is relatively easy to obtain. We only have to notice that 
if p and p' are different multiplets (i.e. they differ in at least one element), then, with obvious 
notation, for every finitely-supported operator 0, 



lim (y(N)\0\V(N!)) = 0. 

N— >oo 



(4.20) 



Indeed, O changes the number of sites with vector Vj (see (2.1)) only by a finite amount, for 



every j. Since for different p and p' , there are j's for which Nj and N'- differ by arbitrary large 



amounts as N — > oo, the result follows. Along with M.10), this yields 



lim Vat = V \c p \ 2 ^(p). 

N—>co ' — ' - — 



(4.21) 



5 Asymptotic behaviour of the entanglement entropy: subman- 
ifold cases 



We now analyze Formula (2.11) and Conjecture [3j first providing verifications for various cases 
where the state ip corresponds to an integration on a submanifold of <CP D , then generalizing the 
analysis to a saddle point argument for general submanifolds. 



5.1 Finite linear combinations 



For finite linear combinations, we have found the measure (4.6) supported on the finite set V 
with weight |c v | 2 - In such states, the Renyi entropy of a block of spins of length m is, according 



to (2.11), 



S n = log 

1 — n 



i 

E 

. v Q GV 
\ q: — 1 , . . . , r< 



n 

a=l 



IK- 



a+1 



a=l 



(5.1) 



(recall that the vectors v a have length one). Like for the spin | case [TT], we find, for all spins, 
saturation at m large 



lim Sn 

m— >-oo 



(5.2) 



That is, the entanglement entropy of the thermodynamic limit of a finite linear combination of 
basic zero entropy states |^ v ) reaches a finite maximum as m — > oo. This maximum is obtained 



from the sum in (5.1) by taking the terms where all vectors v a are aligned. 
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5.2 Elementary states 

We now present a verification of Conjecture [3] in the case of the permutation symmetric states 
(2.22), using the result expressed around ( |2.23| ) and derived in Subsection 4.2 In this and the 
next subsection, we will use the notation 



D 



v (m : = n 



\k=0 



Poe 



Pie 



ie 



with 



D 

52 p* 

j=0 



1. 



(5.3) 



/p], and with an additional overall 6 and p-dependent phase 



This is essentially ^ v (g,a) w ^h a j 
that we put for convenience. 

Thanks to (4.11), we effectively have a (continuous) infinite linear combination of the zero- 
entropy states (4.18). The fact that the linear combination is now infinite allows the entangle- 
ment entropy of a block to diverge as the size of the block increases (rather than saturate, as in 
the previous case) . According to Theorem [2j the entropy in the large- N limit is simply 




n 



d D e a 



IIv t (f,P)-v(f +1 1 p) 



a=l 



(5.4) 



where we have taken the measure ^ given in (2.23) (and the additional phase in (5.3) cancels 
out). Note that, in the scalar product above, the variables p are fixed and identical for every 
copy a. Thus, the particular infinite linear combination of zero-entropy states considered in this 
example only involves the subset of states which is parametrized by the same values of p, that 
is states with fixed magnetization. For such states the scalar product of single site vectors is 
given by 



D 



D 



A(a a 



,p)-v(e a + 1 ,p) = l[e- i W 



a + l_ 



!)p* 



52 pj 



i(6° 

o V 7 



n 



(5.5) 

fe=o 3=0 

As previously done for the s = 1/2 case in [11] , we will now carry out a detailed saddle point 
analysis on the integral (5.4), elaborating on the results presented in |llj . First, we expand the 
logarithm of the scalar product above about = for all j and a 



D 



log n 



k=0 



D 

52 Pi 

j=0 



i(0 



a + l 



\52p*pM +1 

3^k 



+ 



1 



D 



2 ^ 

k=0 



Pk(p k -iW k 



lCt+1 



(5.6) 



neglecting terms of order 3 and above. Note that this is purely real and has no first-order term; 
this is thanks to our choice of phase factor in (5.3). This yields 

n 

H [ v t(r,p)-v(r +1 ,p) 



a=l 



(5-7) 

As we can see from (5.4) and (5.7) although the scalar products depend on the variables 9] 



D 



these variables are not integrated on in ( |5.4| ) (recall the definition (2.6)). They represent the 
fact that the vectors (5.3) are fixed up to a phase. We choose once more Of) = 0. 
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Let us plug (5.7) into (5.4) and change variables to t" 



mpj{ 



9? 



9j) for a = 2, . . . , n and 

tj = Oj for j = 0, . . . , D — 1. With this change of variables the integrand does not depend on the 
variables tj so that D of the integrals can be computed, contributing a factor (2n) D which is 
canceled by the same factor in the denominator. The remaining integrals are of Gaussian type 
and can be computed using standard techniques: 



n 



log 



777 



„ \ 1— n 
D-l \ n 

3=0 J a=2 y ' 



d D t a e # 



k fc=0 K 



xe 



n-l 

E 

a = 2 



(5.8) 



Therefore, 



with 



D-l 



S n ~ S log(4vr 2 m) + log I ft Pj I + 



(5.9) 



j=o 



F{n) = JJ / ^ e 



a=2 



E 



15-1 _ 1 



n-l 

-E 

a=2 



(5-10) 

This expression clearly shows the leading s log m behaviour of the entropy for m large already 
given in [6]. The additional terms can be obtained by computing the Gaussian integrals using 
the general result 



(2vr) M p |BM-i.B 
detA 



(5.11) 



where, in our case, the number of variables is M = D{n — 1) and B = 0. The det A can be 
computed in general to: 



det A = n D (- ^ (D - 1 - k)a k + a D \ 
V fc=o / 



n-l 



(5.12) 



where are elementary symmetric polynomials on the variables { , . . . , ^ x 1 } • h> ue to 
the property J2iLo Pi = 1 it turns out that the expression above is equivalent to, 

/ \ n-l 



det A = n 



D PD 



nD-l 
. 3=0 Pi. 



and therefore, 



F(n) 



n-l 

'IS* pA 2 



(2vr) s ( 1 - n )n s I p D 



(5.13) 



(5.14) 



This immediately gives (2.24) and (2.25). The expression (2.25) exactly agrees with the result 
obtained in [Tj. These results may also be compared to the expressions found for the s = \ 
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case and the particular state of the Ferromagnetic XXX chain considered in |10| . We referred to 
this state as being associated to a "great circle" on the two-dimensional sphere (Bloch sphere). 
In the present language, the state in question is the zero- magnetization permutation symmetric 
state corresponding to setting Po = Pi = \. For these particular values we recover the results 
given in equation (5.73) of [10], as expected. 



5.3 Integration on CP 



D 



The states considered in the previous subsection exhibit logarithmic scaling of both the Renyi 
and von Neumann entanglement entropy as slogm for large blocks of size m. Reviewing the 
methodology employed there, it is straightforward to generalize this result to other linear com- 



binations of the basic states (4.18), with integrations over 9 and summation over p (recall that 

,2 



a i = Pj)- A simple example of such states are the states (4.21) of Subsection 4.3, obtained as 
the thermodynamic limit of a finite linear combination of states (2.1). Another simple example 
is obtained if we choose a number r of the parameters 9j to be fixed. This would reduce the 



number of integrals involved in (5.1 ) from Dn to (D — r)n. Thus, carrying out a similar change 
of variables as in the previous section for the remaining variables the entropy would now diverge 
as (s — |)logm. We did not obtained these states as thermodynamic limits of sequences of 
finite- N states, but the second part of Theorem [T] guarantees that they do occur. Given that 
r = 0, . . . , D this provides a large number of different scaling behaviours which the entropy will 
exhibit, depending on the choice of the initial state. These cases correspond to the measure ji 
in (2.10) being supported on various submanifolds, and not only part I of Conjecture [5] hold for 
these cases, but also Part II, although this requires some more analysis (see the general analysis 
in the next subsection). In addition, if we further give different weights to the various basic 
states |VV(0,p)) in the linear combination, we will find that this does not change the leading 
behaviour of the entropy for large m but will change next-to-leading (constant) corrections and 
re-dependent corrections, in the case of the Renyi entropy. 

More generally, it is also possible to construct linear combinations of the states |VV(0,j>)) whose 
entropy for large m is larger than in any of the cases considered in the previous paragraph. This 
may be achieved by allowing, in the thermodynamic limit, the parameters pj, which where fixed 
in the previous section, to vary in a continuous fashion (i.e. by taking infinite linear combinations, 
or integrals, over the pj as well). We have not shown that these occur as some thermodynamic 
limits, but Theorem [T] guarantees it. If all D = 2s independent parameters are allowed to 
vary one can show that the entropy will then scale as 2s log m. By a similar argument as that 
above, if r parameters pi are fixed the large-block entropy will then be dominated by a term 
(.D — !) logrre. We can therefore see that playing with the number of parameters 9 and p which 
are allowed to vary, the entropy scales as | log m with d = 0, ... ,4s. 

We now generalize the analysis of the previous section to the situation when both sets 9 and 
p are integrated over. In addition, we will consider the most general situation, when we have an 
arbitrary linear combination of states, parametrized by functions g(9 a ,(j) a ) > 0: 



5,1 



1 



log n 



\a=l 



2tt 



d D 9 a 



n 

0=1 



vH9 a ,6 a )-v(9 a+1 ,6 a+1 ) 



(5.15) 



where S D is the D-dimensional sphere. Here, instead of using our original variables p a we have 
chosen to use a natural parametrization of the sphere in terms of angles cp a as given in (2.17) 



with the restriction (2.18) and up to the usual identification p°- 
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Comparing (5.15) to the general formula (2.11), we have that the measure \x has been chosen 
so that 

dfJt(v@., $) = 9(1, 0) d D 0d D ^. (5.16) 
Our starting point as before is the scalar product of single normalized vectors 



D 



D 



)-v( 



.a- 



x )=n 



k=0 j=0 



(5.17) 



v 

Expanding the logarithm of this scalar product about 0? = 0J +1 and </)j = (f>J +1 for all j and a 
we find, up to second order, 

ior [ V t(r,<n 



sin 



2 (1 I (±a+l ia\2 



Y-iuj a (5.18) 



fc=0 



fc=0 \p=0 



where we used 



E 

i=o 



\2 D-l fk-1 



4p? 



En 

fc=0 \p=0 



sm 



vr 1 -^) 2 , 



(5.19) 



valid to second order. In this expression we have used that to second order, the coefficients of 
the quadratic terms above can be written in terms of the variables only (or equivalently 4>j ) , 
for each value of j. The real quantity co a is given by 



D 

E 

3=0 



1 



(5.20) 



Contrary to the case of the previous subsection, the pure imaginary part — iu a cannot in 
general be canceled by a choice of phases of the vectors v(0 a , <p a ). In particular, the quantity 



U) ' 



(5.21) 



a=l 



is not zero for any n > 3 (it is zero for n = 2). However, since the first term under the summation 
symbol in (5.20), when summed over a, is invariant under 0® \— > 0® + qj for any qj, the quantity 
is in fact a function of 0" +1 — 0^ and +1 — (for j = 0, . . . , D, a = 1, . . . , n) only. From 
(5.20) it is clear that it is quadratic in these differences. 



Putting everything together we have 

n 



aQ ^ a )-v(0 a+1 ,6 a+1 ) 



E 

ce=l 



j^k fc=0 



a=l 



n -D-l / k — 1 

e e n « 

a=l fc=o \p=0 



-f E E I n sin 2 ^ (^ +1 -^) 2 



(5.22) 



We may now substitute this expression into (5.15). By symmetry under exchange of copies of 
the integration measure in (5.15), we find that the result of the integral in (5.15) is real. Hence, 
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we may replace, in the integrand, the factor e min by its real part cos(mO) ~ e m n / 2 (in the 
saddle point approximation). Since £1 is quadratic in the differences — Of and p" +1 — p'j, 
then Q, 2 is quartic. This shows that the phase factor e~ mt ^ in (5.22) can be neglected in the 
saddle point approximation of the integral. 

We now proceed as in the previous section. The integrals in the variables 9® may in fact be 
carried out exactly as before, with the only difference that the parameters are not constants 
in the present case and therefore their product must be kept inside the integral. After the change 



of variables = ^Jmp x -{9°- — 6>j) and 
have: 



b}) for a > 2 and t) 



9 1 h 1 



we 



2 s log 



m slogn 1 

-= — H log 

2tt 1 — n 1 — n 



l-n 
2 



n-1 D-l /fc-1 



\[dh c 

a=2 



a-2 fc=0 \p=0 



/ V / dWlIp} aii^T (5.23) 

JO J S D \ ,_n / 

\ 



D-l I fc-1 



-|E E H sin 2 ftl (^ +1 -^) 2 -| £ n sin 2 ^ ((/x 2 J 2 +(^) 2 ) 



fc=0 \p=0 



/ 



The integrals in for a > 2 can be easily carried out as they are of the standard Gaussian 
type. Moreover, there is factorization for each fixed j. Integrating we obtain: 



(27T 



,3(71-1) 



11 



D-l fc-1 

n n 

k=0 p=0 



n-1 • 



(5.24) 



n nsin^ 



substituting into (5.23) after expressing the product n^Lo^) ^ n terms of the angles hj we find 
m 2s log n 



2s log 



+ 




(D-l fc-1 

dtfg&h 1 )" H ^sin(2/4)n sin2/l : 



1— n\ 



2 ul 
V 



(5.25) 



fc=0 



p=0 



Hence, the von Neumann entanglement entropy takes the form 



Sl = 2slog^- dt 1 [ dh l g(t\h l )\og{ . . - 9{ -\- \ , 



sm 



2 K 



(5.26) 



where we have used the fact that 



r2n r 

/ dt 1 / dh}g(£,h}) = 1. 
Jo Js D 



(5.27) 



In order to compare this result to the general formula (2.14) we need to find the relationship 
between our function g{t^,h}) and the function defined in Conjecture [3J1I. Recall that the 
defining property for this function is 



dn(v(t,h)) = d 2D vf(t,h) 



(5.28) 
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and d 2D v is defined in (2.20). Therefore, 

f(t\h 1 ) = - 



g(t\h U 



and requiring that 



2^ D n 1 sin(2/4)n 1 sin 2 h l 

k=o p=0 



(5.29) 



/f'Z7T /* 
d 2D vf{t l ,h l ) = dt 1 dh}g(i,h x ) = l, 
Jo Js D 

the von Neumann entropy becomes 

Si = 2* log |^ - J ^v/^^log/U 1 ,^ 1 ), 



(5.30) 



(5.31) 



as expected. Similarly, the Renyi entropy is that obtained from (2.13) with d = 4s. 



5.4 General case 

The argument in the general submanifold case is a consequence of a similar saddle-point approxi- 
mation as that done above. Let Xj be real coordinates on the submanifold (formally representing 
all coordinate patches), with j = 1, . . . ,d where d is the dimension of the manifold. Then the 
state that we are considering is 



tp = d d x \J r](x) f(x) ip 



v(x) 



(5.32) 



where r)(x) is the determinant of the metric, at the point x on the submanifold, induced by 
the Study- Fubini metric on CP D ; and where v(x) is the homogeneous coordinate on CP^* 



corresponding to the point x. Hence according to (2.11) we wish to evaluate, in the large-m 
limit, the quantity 



Sn 



1 



1 — n 



log 



f (fld d X«y/^f(x<*)) I f[ vt(x Q ) • v(s 
J \a=l J \a=l 



(5.33) 



Note first that (5.22) implies, to quadratic order (with assumed normalization v] • v c 



a=l 



a=l 



-in 



(5.34) 



where £1 is quadratic in the difference of the coordinates of v a and v Q+ i on CP D (these are 
coordinate-independent statements). Further, the fact that the integration measure in ( |2.11 ) is 
symmetric under exchange of copies implies, as in the argument after (5.22), that e™^ can be 



neglected under the integral in (2.11) in the saddle point approximation. Hence we may replace 



in (2.11) 



n 

a=l 



(5.35) 



where we use the distance on CP^ given by (2.16). The integration in (5.33), in the large-m 



limit, is then performed along entirely similar lines as those of the previous subsections. We use 



4) 
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(to quadratic order), where rj^ k {x}) is the metric on the submanifold, here evaluated at the 



coordinates x 
a > 2 and ij ■■ 



belonging to the copy a = 1. After the change of variable t" = ^Jmi^x - 
x], the result of the Gaussian integration over if, a > 2 is (2.13). 



6 Asymptotic behaviour of the entanglement entropy: fractal 
subsets 

In this section we complete the verification of Part I of Conjecture [3] by providing an explicit 
example of a case where ip corresponds to a measure over a fractal subset of CP^ with D = 1 
(i.e. spin s = 1/2). We then give the general argument for integrations over fractals that was 
briefly presented (in the special case s = 1/2) in |11| . 



6.1 A Cantor set on the Bloch Sphere 

Let us consider the entanglement entropy for s = \ and provide an explicit example in which the 
geometric dimension d is a fractal dimension. In our example this will be the fractal dimension 
of the well-known Cantor set (see e.g. [29]). Here we identify the Boch sphere S 2 with CP 1 in 
the usual way, via 



z e 



i0 



x + iy 



z 2 e ie 



where (x,y,z) £ S 2 C R 3 , and v is a homogeneous coordinate on CP 1 . 



Consider once more our expression for the entropy (2.11). Let |c v | > for every v in a 



finite subset of unit vectors, and consider the associated quantum state 



(6.1) 



Its entanglement entropy is 




n 



(6.2) 



where, as in the previous subsection, we used the normalization |v a | 2 = 1. Note that the 
argument leading to (5.35) holds true for any measure of integration. Hence, in general to 
quadratic order, we may make the replacement (here specialized to s = 1/2) 



n 

a=l 



v+ • v Q+ i !-)• exp 



1 



n 
a=l 



V a +1 



(6.3) 



where v a G S 2 is the coordinate on S 2 representing the homogeneous coordinate v Q G CP 1 . 

Consider instead a state in a linear combination supported on the Cantor set (the method 
considered here can be applied to other similar fractal sets). For defmiteness, we consider a 
Cantor set on half of a great circle on the Bloch sphere, as shown in Fig. 1. The great circle 
that we consider is that with z = 0. The Cantor set C on the unit interval [0, 1] is the limit 



C 



lim Ck 

k— >oo 
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Figure 1 : A Cantor set (pointed by the blue arrows) on half of the equator (green circle) of the 
Bloch sphere (semi-transparent red) . The points are computed by employing the usual iterative 
definition of the Cantor set up to 5 iterations, and the map (6.4). 



where the finite sets Cj- : k = 0, 1, 2, . . . are defined recursively by 

C k+ i = * (C fc U (2 + C k )) , Co = {0} 

(here we use the general notation a + b{q\, q2, . . .} = {a + bqi, a + bq2, ■ ■ •})• We embed it on the 
Bloch sphere via the map <fi : [0, 1] — > S 2 , defined by 

4>(x) = (cosnx, sin7rx, 0). (6-4) 

The quantum state V'cantor i s a Hausdorff integral over the set 0(C) C S 2 ; for simplicity, we take 
it with constant coefficients. We denote the Hausdorff measure on C as /j,q. Hence, we consider 

V'cantor = / c ^/ i c( v ) with measure normalized as / dfi c (v) = 1. (6.5) 

J<t>{C) J<j>{C) 

The measure \i c can be defined recursively following the recursive definition of the Cantor 
set itself. We note that in the finite set C k , at step k of the recursive definition, there are exactly 
2 k points. Hence the step-fe quantum state ?/> fc associated to C k with constant coefficients is 

= ( 6 - 6 ) 

ve^(C fc ) 

Then 'i/'Cantor is defined by taking the limit: 

V'Cantor = , lim (6.7) 
fc— >oo 

One can show that this limit exists weakly (i.e. in the local-operator topology on infinite-iV 
(thermodynamic) states). 

Consider now the entanglement entropy S n = S n (ijj) as a function of the state ip. It is clearly 
continuous as it is a continuous function of the evaluation of -0 on a finitely-supported operator, 
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see (3.3). Hence in order to evaluate the entanglement entropy of V'cantor) we mav use formula 
(5.1) and limS n (ip) = S n (limip), giving 



Cantor) ~ , nm S n (lj} k ) 
k— »oo 



1 — 71 



log I lim 

k— >oo 2 



- y 

kn / j 



{v a &<t>(C k )} La=l 



(6.8) 



At large m, the only terms that give important contributions are those for which all overlaps 
IVv a+:L ) (f° r an a ) are near to 1. When such an overlap is near to 1, we can use the 
approximation (6.3). Further, for two vectors v = <j)(x) and w = 4>{y) near to each other on S 2 
(i.e. with \x — y\ small), we can approximate \v — w\ 2 ~ ir 2 \x — y\ 2 . Hence, for the purpose of 
evaluating the large-m leading logarithmic behaviour of 5n (V'Cantor) we can use 



E 

{v Q G0(C fe )} 



.0=1 



£ ex p 

{x a £C k } 



1 n 

it m 



a=l 



(6.9) 



Now we simply use the recursive definition of Ck and the fact that the large-/c limit exists. 
We have 



lim 7T" 



exp 



9 ?1 

7T m 



{x a £C k } 

{x a ec fe+ i} 



E 



9 n 
TT m 



E 

0=1 



— lim — j- 
2 n fc^oo 2 fcn 



E 



exp 



{x a &±(C k U(2+C k ))} 



9 n 
Tim 



E 



1 1 

— urn — j— 
2 n fc^oo 2 kn 



E 



exp 



vr 2 (m/9) 



a=l 
n 



1 - x c 



a=l 



3.10) 



{x Q ec fc u(2+c* fc )} 

For the large-m power law asymptotic, we can write in ( |6.10[ ) 

E « E + E • 

{x a <=C k U(2+C k )} {x a eC k } {x a £2+C k } 

Indeed, whenever a nonempty proper subset of {x\, . . . ,x n } is in Ck and the rest is in 2 + Ck, 
there will be finite contributions in the exponential which make it decay exponentially at large 
m. Since both sums then give the same contribution, we obtain 



Qr, 



1 



lim 



2"-i fc^oo 2 k 



1 



i £ ex p 



7T 2 (m/9) ^ 2 



a=l 



3.11) 



Defining qi = log 2 n-i Qge (for £ large) we then find the recursion relation 

qi+l = qi-1. 



3.12) 
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Since this must be true for a continuum of £ (because we are looking at the asymptotic region 
where the integer 9^ is large) and since qg is monotonously decreasing with £, the unique solution 
is 

q e = r-£ (6.13) 

for some number r. Hence, 

_ n-l log 2 

Q m = Rm 2 (6.14) 

(for m large) for some other number R. 
From this, we obtain at large m 

all \ 1 1 ( -^^A llo g 2 ! ^Cantor, ,„,,, 
SWlV'Cantor) ~ Y~n g V™ ) = 2 fog~3 g m = 2 g ™ ^ ' 



where dcantor = log2/log3 is the fractal dimension of the Cantor set. This agrees with (2.12). 



6.2 General argument 

As is clear form the calculation above, a crucial feature of the Hausdorff measure that is used 



in obtaining (2.12) is its scaling property. 

Let v be some probability measure on M. 2D which is absolutely continuous with respect to 
the d-dimensional Hausdorff measure 7~L d on M. 2D . Hence, it is supported on a fractal set of 
dimension d, and can be written 

dv(x) = f(x)dH d {x) (6.16) 
where x G ]R. 2D . Note that the Hausdorff measure satisfies 

dU d (sx + y) = s d dU d (x) (6.17) 

for all x, y G M. 2D , s > 0. Let us assume that lim s _;.o sup(/(x) : \x\ < s) = M is finite. This 
implies that, if F is a positive function on M. 2D that decays fast enough at infinity, then 

< J dv{x) F(x/s) = s d j dU d {x) f(sx) F(x) < Ks d (6.18) 

as s — > 0, where K = M J dH d (x) F(x) is finite. 

Let v G CP D (homo geneous coordinates), and let U v C M? D be a coordinate patch around 
v (which can be see as an open subset of the tangent plane at v). Let <p v : U v — > CP D be the 
associated coordinate diffeomorphism, with </> v (0) = v. Given any w G (j) v (U v ), we can define 
x = (/>~ 1 (w) G U v . We normalize 4> v in such a way that if w and w both tend to v, then 

D(w, w) ~ \x — x\. (6.19) 

Let \i be a probability measure on CP D such that for every w G CP D in a neighborhood of 
v, there is a measure u v as above such that 

d/j,(w) = di/v^-^w)) = du v (x). (6.20) 

This is a measure supported on a fractal set on CP^ (of fractal dimension d) , and by integration 
over a fractal set on CP^, we mean integration over such a measure. 
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Using (2.11), (5.35) (again noting that this holds for every measure \x in (2.10)), (6.19) and 
(6.20), we have 



n 



log 



1 



n 



log 



[(nM**)) (f[e 

J \a=l / \a=l 

n 

n dv vi (. 



-f D(v Q ,v Q+ i) 




a=2 




n-l 

J2 \x a + 1 -x a \ 2 +\x 2 \ 2 +\x n \ 2 

a=2 



a=2 



Since fi is a probability measure (hence has a finite integration), it must be that for almost 
all v with respect to fj,, the measure dv v (u) has an associated density f(u) that is finite in a 
neighbourhood of u = 0. Hence, we can use (6.18) (multiple integral generalization) to find 



< 




n-l 

J2 \x a+1 -x a \ 2 +\x 2 \ 2 +\x n \ 2 

a=2 



d(l-n) 

< K (vi) m 2 



a=2 



where if(vi) is almost surely finite. In order to arrive at the claimed asymptotic behaviour, we 
must guarantee that there is a subset of /i-measure nonzero for vi such that the lower bound 

d(l— n) 

for the a bove integral is strictly greater than 0, and further of the order m 2 . In the context 
of (6.18), this occurs if the limit lim s _ > o/(sx) = f(x) is integrable and integrates to a nonzero 
value. It seems natural to believe that this holds under i/ v -integration for almost all v with 
respect to fj,; it in fact holds for densities f(x) that are continuous on their support. Assuming 
this lower bound property, we then find 



S n ~ ^ log m 



3.21) 



in agreement with (2.12). 



7 Conclusions and outlook 

In this paper we have studied the von Neumann and Renyi bi-partite entanglement entropies 
in the thermodynamic limit for spin-s quantum systems. This includes quantum spin chains 
as well as any other quantum many-body systems whose Hilbert space is described by a tensor 
product of spin-s sites. In our study, we have considered a very special type of quantum states 
known as permutation symmetric states. These are states which are permutation symmetric 
under the exchange of any sites. In the present context, this type of states was considered in [7] 
and part of our work has been devoted to extending the results of this work and to set them in 
a more general context. 

The motivation to study the entropy of permutation symmetric states is two-fold: on the one 
hand this kind of states, for finite chains, are closely related to Dicke states, which are widely 
studied in quantum computation and have been experimentally realized [30]. On the other hand, 



the behaviour of the entropy for large subsets A, as defined in (1.2), is unusual and not dictated 
by conformal invariance. As for conformal critical points, logarithmic scaling is found. However, 
unlike in the latter theories, the numerical coefficient of the logarithm is not related to a central 
charge and is the same for both the von Neumann and Renyi entropies (e.g. n-independent). 
Hence in particular, it is also the same for the so-called single-copy entropy (the limit n — > 00 
of the Renyi entropy). In this paper, following on previous work [TQl [Tl] for the spin-| case, we 
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have shown that the coefficient of the leading logarithmic behaviour has a geometric meaning; 
it measures the geometric dimension of the set of basic zero-entropy states in terms of which the 
permutation symmetric states can be expressed. For spin-^ the geometries involved are subsets 
of the Bloch sphere. For spin-s we have shown that this generalises to CP 2s . 

From a mathematical standpoint the most challenging aspect of our work has been to provide 
a meaningful prescription which would allow us to carry out the thermodynamic limit, that is 
when the length of the chain iV and the dimension of the Hilbert space tend to infinity. We have 
provided plausible arguments as well as examples tending to prove the existence of a probability 
measure \x supported on all subsets of CP 2s with respect to the Study- Fubini metric (the natural 
metric in CP 2s ). Provided this measure exists (Theorem 1), we have proven that correlation 
functions of the infinite chain as well as its von Neumann and Renyi entropies may always be 
expressed as integrals over such measure in a very general form (Theorem 2). 

When considering the thermodynamic limit of permutation symmetric states we have in 
particular considered a very particular basis for such states, namely that previously introduced 
in the work of Popkov et al. [7] as a generalization of their previous work [6] for the spin- 
I case. In this basis, which we referred to as the elementary vectors or states, permutation 
symmetric states are characterized by fixed ratios pi of spins with S^-projection given by s — i 
and i = 0, 1, . . . , D. In [7] the logarithmic scaling behaviour of the von Neumann entropy for 
these states had already been found. Here we have reproduced and extended that result to 
the Renyi entropy, giving the scaling behaviours (2.24) and (2.25). Popkov et al. were able to 
obtain their results by exploiting the particularly nice combinatoric features of their chosen state 
such as the fact that non vanishing eigenvalues of the density matrix are generalised binomial 
coefficients. For large blocks, a saddle point analysis of these coefficients then leads to the 
behaviour (2.25). In our case this behaviour follows from Theorem 2, which gives our particular 
formula for the entanglement entropy based on twist operators [10]. Further, the elementary 
vectors are very particular linear combinations of our zero-entropy vectors, and thanks to the 
generality of this formula we have been able to investigate the entanglement entropy of much 
more general linear combinations (both finite and infinite) and we have shown that the scaling 
behaviour of the entropy will generally be different depending on the chosen state. Besides the 
case above, we have explicitly computed this behaviour for the state described in Subsection 
5.3, giving (5.31). In this case, the coefficient of the logarithmic term is 2s and is maximal, 
meaning that no linear combination of permutation symmetric states can have a faster growing 
entropy than this particular one. We have also explained how other states may be systematically 
constructed giving coefficients ^ = 0, |, 1, |, . . . , 2s for the logarithmic divergency, and how in 
these cases the O(l) correction term to this divergency has an interpretation in the context of 
geometric quantum mechanics. Finally, we have shown how one may construct states for which 
the quantity ^ takes irrational values, with d representing the fractal dimension of some fractal 
subset of CP 2s . We have provided an explicit construction for s = 3 and the Cantor set and a 
general argument for higher spins and general fractal sets. 

We leave to future work a more extensive discussion and proof of Theorem 1, that is the 
existence of a measure /i with the given properties. A rigorous proof of this statement turns out 
to be quite involved and is closely related to the famous moment problem in mathematics which 
asks the question: given a set of moments xq,xi,X2, ■ ■■ under which conditions is it possible 
to find a measure /U such that x n = J x n dfi? A general answer to this question in the higher- 
dimensional case is not known; however in our particular application we believe that such answer 
can be found. 

It would be interesting to investigate further the precise relationship between our results and 
the geometric entropy defined in the context of geometric quantum mechanics. Also, a numerical 
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test of our findings would be of interest. 
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